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ABSTRACT 

In this Letter we investigate the role of recombination energy during a common envelope 
event. We confirm that taking this energy into account helps to avoid the formation of the 
circumbinary envelope commonly found in previous studies. For the first time, we can model 
a complete common envelope event, with a clean compact double white dwarf binary system 
formed at the end. The resulting binary orbit is almost perfectly circular. In addition to con¬ 
sidering recombination energy, we also show that between 1/4 and 1/2 of the released orbital 
energy is taken away by the ejected material. We apply this new method to the case of the 
double-white dwarf system WD 1101+364, and we find that the progenitor system at the start 
of the common envelope event consisted of a ~ 1.5M@ red giant star in a ~ 30 day orbit with 
a white dwarf companion. 
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1 INTRODUCTION 


The formation of a compact binary system composed of two white 
dwarfs (WDs) is widely accepted to include a common envelope 
event (CEE), at least during the last episode of mass exchange be¬ 
tween the first-formed WD and a low-mass red giant (RG). Low- 
mass RGs have a well-defined relation between their core masses 
and radii, providing for DWDs the best-defined state of a progeni¬ 
tor binary system at the onset of the CEE among all known types of 
post-common envelope (CE) systems. As a result, DWD systems 
have served extensively as test-sites for attempts to understand the 
physics of CEEs, using both population synthesis approaches and 
hydrodynamical methods. 

Previous attempts to model a CEE between a low-mass RG 
and a WD did not succeed to eject the entire CE during three- 
dimension al (3D) hydrodynamical simulations (f or most recent 
studies, see lPassv et al.l20l JlRicker & Taaml2012l) . The final state 
of these simulations is that a significant fraction of the expanded en¬ 
velope remains bound to the formed binary, forming a so-called cir¬ 
cumbinary envelope. Then almost no energy transfer can take place 
from the binary orbit to that circumbinary envelope. Observation- 
ally, no circumbinary disk in a post-common envelope system has 
been found so far. 

It has been proposed long ago that recombinatio n energy of 
hydrogen and he l ium should play a role during a CEE i Lucvll967 

I- 


; IPaczvriski & Ziolkowskil Il96a : Han et all 11994 


20021) . However, until now, this energy was not yet taken into ac¬ 
count in 3D modeling. While the initially available recombina¬ 
tion energy can be easily comparable to the binding energy of the 
remaining bound envelope (e.g. [Passv et al.ll2012l) . the important 
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question is when and where the energy is released - to be use¬ 
ful, recombination energy should not be released too early in the 
CEE nor in material already ejected, but instead in the circumbi¬ 
nary envelope at a time when the recombination energy is compa¬ 
rable to the binding energy of the not-yet ejected material. In this 
Letter , we investigate if the inclusion of recombination energy can 
help to remove the circumbinary envelope. We apply the new ap¬ 
proach to the system WD 1101+364, a well-measured DWD that 
has P or b = 0.145 d and a mass ratio of q = M\jMi = 0.87±0.03, 
where M\ ~ 0.3lMg and M 2 — 0.36M, ; , ar e the masses of the 
younger and older WDs, respectively (M ars 319951) . 


2 INITIAL SET UP AND DEFINITIONS 

We anticipate that the progenitor of WD 1101+364 was a low- 
mass RG that had a degenerate core of 0.31A/q. We consider the 
range of masses for the RG donor, j, from l.OM 0 to 1.8 M@. To 
evolve the RG and find the initial one-dimensional (ID) stellar pro- 
file, w e use TWIN/S tar stellar code (recent updates described in 
iGlebbeek et al.ll2003) . The stars are evolved until their degenerate 

He cores have grown close to 0.31M 0 . _ 

For 3D simulations, we use StarSmasher dGaburov et al.l 
l20ld : iLombardf et al .1 [ 20 TTI 1 . a Smoothed Particle Hydrodynam¬ 
ics (SPH) code. Technical details o n usin g this code to treat CE 
events can be found in lNandez et al .1 J20l4) . A ID stellar profile is 
imported to StarSmasher, where an initial stellar model repre¬ 
sented by a certain number of particles N p is generated via a re¬ 
laxation process. The core of a RG is modeled as a point mass - a 
special particle in SPH that interacts only gravitationally with other 
particles. Because the centre of the giant is not fully resolved, the 
core mass, M c \, is slightly more than in the ID code (see Table 
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Table 1. Initial conditions 


Model 


^env,l 

Mc,i 

Rlof 

^orb,ini 

^orb,ini 

Np 

n 

^bind 

E ac 

^orb,ini 

^tot,ini 

A 

1.0RG0N 

0.985 

0.668 

0.317 

28.21 

60.11 

46.57 

99955 

0.00 

-8.454 

2.059 

-1.118 

-7.513 

1.047 

1.0RG1N 

0.985 

0.668 

0.317 

29.78 

63.51 

50.57 

99955 

1.00 

-8.454 

2.059 

-1.058 

-7.453 

0.992 

1.0RG2N 

0.985 

0.668 

0.317 

31.35 

66.81 

54.57 

99955 

2.00 

-8.454 

2.059 

-1.006 

-7.401 

0.942 

1.2RG2N 

1.195 

0.877 

0.318 

29.48 

60.74 

44.00 

99955 

2.00 

-12.328 

2.725 

-1.345 

-10.945 

1.093 

1.4RG2N 

1.397 

1.079 

0.319 

27.73 

55.59 

36.24 

99955 

2.00 

-16.947 

3.369 

-1.715 

-15.293 

1.217 

1.5RG2N 

1.498 

1.179 

0.319 

25.66 

50.82 

30.81 

99955 

1.65 

-20.636 

4.038 

-2.015 

-18.609 

1.267 

1.5RG2NP 

1.498 

1.179 

0.318 

25.66 

50.82 

30.81 

200221 

1.65 

-20.345 

3.697 

-2.011 

-18.659 

1.285 

1.6RG2N 

1.598 

1.275 

0.323 

25.80 

50.54 

29.76 

99955 

2.00 

-22.837 

4.009 

-2.157 

-20.985 

1.312 

1ARGOS 

1.598 

1.275 

0.323 

31.25 

48.61 

27.97 

99955 

0.00 

-22.358 

3.997 

-2.241 

-20.602 

1.106 

1.7RG2N 

1.699 

1.376 

0.323 

22.83 

44.25 

23.78 

99955 

2.00 

-28.638 

4.338 

-2.619 

-26.918 

1.356 

1.7RG0S 

1.699 

1.376 

0.323 

27.58 

42.97 

22.67 

99955 

0.00 

-28.003 

4.326 

-2.694 

-26.371 

1.148 

1.8RG2N 

1.799 

1.481 

0.318 

16.34 

31.37 

13.86 

99955 

2.00 

-44.167 

4.676 

-3.912 

-43.404 

1.401 


The models names are composed as following: two digits representing the RG mass are followed by “RG”, /] value is is outermost smoothing length; 
then “N” stands for non-synchronized and “S” for synchronized (“S”) cases. “P” denotes the model with about twice larger number of particles than 
in all the other models. j, A/ env ,j and M c i are the total, envelope and core mass of the RG, in Af 0 . R,iof' s the radius of the donor Roche lobe, in 
/? . , and 7) describes the adopted donor’s radius definition (see §2). a or bmi is the initial orbital separation in R 0 , P 0 rb,ini is the initial orbital period in 
days. Np is the total number of SPH particles that represent the RG. £bind> E rec , £ or b,ini and P| OIJm are the binding energy of the RG envelope without 
recombination energy, the total recombination energy of the RG envelope, initial orbital energy, and initial total energy , defined as the sum of the 
binding, recombination, and initial orbital energies, respectively, in the units of 10 46 erg. A is a star structure parameter (see Equation[5J. 


Q] for this and other initial values). This ensures a proper match¬ 
ing of stellar profiles of 3D envelopes with ID stellar profiles. The 
envelope mass in a 3D star is M env .i = \ — M c .i. 

In a 3D star, the radius of the star, R$ph, can not be defined as 
uniquely as the p hotos pheric radius of the ID star (for a thorough 
discussion, see iNandez et alJl2014l) . The stellar radius can be pa¬ 
rameterized as Psph = Rom. +1] (touti where R out is the position of 
the outermost particle and /z out is the smoothing length of that parti¬ 
cle. The parameter r] can range between 0 (in this case, some mass 
will be found above Rsph) to 2 (with all mass contained within 
I?SPh). In addition, we note that a synchronized giant is expected to 
attain a larger radius after relaxation than a non-synchronized giant. 

The initial orbital separation, fl 0 rb,inb for the non-synchronized 
cases, is found from the assumption that Rspu is equal to the Roche 
lobe (RL) ove rflow radius, R r i 0 f, and using the approximation by 
lEggletonl Jl983h . The initial orbital period, Rorb.ini is found assum¬ 
ing a Keplerian orbit. For the synchronized cases, the orbital period 
and separation are founded the moment when the RG overflows its 
Roche lobe (see §2.3 of lLombardi et all201 lh . 

Equations of state (EOS). The standard EOS (SEOS) in 
StarSmasher is analytical and includes radiation pressure and 
ideal gas contributions. To take into account recombination energy, 
we need another prescription for the EOS. Because we evolve the 
specific internal energy Uj and density p, for each particle (among 
other variables), we prefer an EOS that uses m; and p; as indepen¬ 
dent variables. However, such an analytical expression does not ex¬ 
ist in simple form when we consider recombination/ionization of 
atoms. Therefore, we are bound to use a tabulated EOS (TEOS) 
which uses u ; and p, to provide the gas pressure Pgas,i, temperature 
T{, specific entropy s etc. 

We use the MESA-EOS module to calculate such tables (see 
§4.2 of lPaxton et al .11201 lh . The core of a RG is modeled as a point 
mass, and the rest of the star has uniform composition. Hence, only 
one table with a single set of composition for H, He, and metals 
needs to be generated for each RG. The tables that we generate 
operate in 9.84 ^ log«[erg g 1 ] ^ 19.0 and —14 ^ logp[g cm~ 3 ] ^ 
3.8. When a particle has a density or specific internal energy outside 
the limits of our tables, we switch to the SEOS. 


Energy formalism. The energy formalism compares the 
donor’s envelope binding energy Ebind with the or bital energy 
before, E nr h.ini, an d after the CEE, E 01 -b,fi n dWebbinkl 119841 : 
iLivio & Sokedl 1988tl : 


(-bind Q-bind((t 0 rb,fin (-orb.ini) — O-bindA^orb • (1) 


Here 0!bind is the fraction of the released orbital energy used to 
expel the CE, 0 ^ Ofbmd ^ 1- The binding energy of the donor’s 
envelope, in its standard definition, is 


(-bind — ^ m i 



3 kTj 
2 p, w H 



( 2 ) 


where m;, 7), p,-, p, and fl ; are the mass, temperature, density, spe¬ 
cific gravitational potential energy, and mean molecular mass, re¬ 
spectively, for each particle i. The constants k, a, and mu are the 
Boltzmann constant, radiation constant, and h ydrog en atom mass , 
respectively, while p; is calculated as in lHemquist & Katzl (1 1989h . 
Note that Ebind in its standard definition does not include recombi¬ 
nation energy. 

The binding energy of the donor’s enve lope is frequen tly pa¬ 
rameterized using a parameter A, defined as lide Koollll990lf 

^ tft'id. l -Hein | 

(’-bind^rloi 

Here G is the gravitational constant. For low-mass giants, A is 
known to have a value close to one, and we obtain similar results. 
We find the orbital energy of the binary system according to 

£orb = lv\Vn\ 2 + Y t ^m i <l>i-Y t ^mj<l>J Li (4) 

i j k 


where p = M\Mi!(M\ -VMi) is the reduced mass, and V 12 = V\ — 
V 2 is the relative velocity of the two stars. The first term gives the 
orbital kinetic energy. The second term is the total gravitational 
energy of the binary, with the sum being over all particles l in the 
binary. The third and fourth terms correspond to the removal of the 
self-gravitational energy of the donor (the sum being over particles 
j in star 1) and of the WD (the sum being over particles k in star 
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2, initially only the WD), respectively: the remaining gravitational 
energy is then just the orbital contribution. 

Recombination energy. In our treatment, the internal energy 
provided by the TEOS includes contributions from ideal gas, radia¬ 
tion, and the recom b inatio n energy for H, He, C, O, N, Ne, and Mg 
(see §4.2 of lPaxton et aljEzoIlh . The recombination energy can be 
extracted as 

„ / 3 kTj aT 4 \ 

Erec — / I ttj — 1 = CTree AE or b, (5) 

i \ 2 Ri' w H Pi J 

where u, is the SPH specific internal energy of particle i. Values of 
£ rec , as expected, scale well with the mass of the envelope. Note 
that here we introduce important new parameter, CT rec - the ratio 
between the recombination energy and the released orbital energy. 

Total energy. The initial total available energy, £tot,ini> is 

£tot,ini = ^orb.ini + ^bind + ^rec • (6) 

This quantity is conserved during the evolution of all our models. 

Bound and unbound material. For each particle, its total 
energy is defined as £ t ot,/ = 0.5m,-v? + m ,</>,- + nijUj, where the first, 
second and third terms are the kinetic, potential, and internal en¬ 
ergies, respectively. If the particle has negative total energy, it is 
bound to the binary. In this case, if the particle is located outside of 
either RL, the particle is in the circumbinary region. Accordingly, 
we classify the particles to be in (i) the ejecta, m un b - the particles 
with positive energy, (ii) the circumbinary, m c ^ - the matter bound 
to the binary but outside of the two RLs. and (iii) binary, - the 
particles are inside either of the two RLs. 

The total energy of the unbound material at infinity is com¬ 
puted when the unbound mass is in a steady state after the CEE: 

^tot.unb = = — a unb^orb • <7) 

i 

Here we introduce oQ b - the energy taken away by the unbound 
material in units of the released orbital energy. Note that in the 
standard energy formalism this quantity is always assumed to be 
zero. 

Angular momentum budget. We calculate the orbital angu¬ 
lar momentum 

•4>rb = 1^12 X ^12, (8) 

where R 12 = Rl — Rl is the displacement from star 2 to star 1. We 
note that the magnitude ,/ 01 -b = |./ 01 -b, z |, where the z-direction is per¬ 
pendicular to the orbital plane. An outcome of a CEE can be charac¬ 
terized by how much orbital angular momentum is lost. We prov ide 
the y-parameter dNelemans et al.l200(il : lNelemans & Tonti'200 5i as 
a way of quantifying angular momentum loss in our simulations: 

M\ + A /2 -/orb.ini -forb.lin ... 

r= ~M - i -• (9) 

^unb •'orb,ini 

3 FORMATION OF A DWD THROUGH A CEE 

Comparison between the two EOSs. We compare the results of 
simulations with our two different EOSs using the model 1.5RG2N. 
In both cases, the initial relaxed stars have SPH profiles that match 
very well the ID stellar profiles for pressure, density, and gravita¬ 
tional potential. However, this is not the case for the specific inter¬ 
nal energy u (see FigureQJ: clearly only the TEOS model matches 
the desired ID stellar profile. As expected, the mismatch between 
the relaxed profile with SEOS and the stellar one is indeed due to 
neglecting recombination energy. 



m [M 0 ] 

Figure 1. Specific internal energy u profiles for the model 1.5RG2N. The 
black asterisks and gray triangles correspond to relaxed 11 profiles for TEOS 
and SEOS, respectively. The black solid line corresponds to the u profile 
from the stellar code. 

We find that the SEOS fails to unbind the envelope in our CE 
simulations. Only about 50% of the envelope becomes unbound: 
the circumbinary matter does not interact with the formed binary 
at all, making it impossible to eject the entire envelope. This re¬ 
sult i s consistent with the fi ndings of previous studies (IPassv et al.1 
l2012l : lRicker & Taamll2012h . On the other hand, the TEOS simula¬ 
tion clearly makes use of the recombination energy and ejects the 
envelope entirely. For all other models presented in the Letter, we 
use the TEOS. 

Masses. At the end of the simulations, we form a binary con¬ 
sisting of M\ and M 2 (see Table [2] for all the outcomes). The un¬ 
bound material M un b is at least 99.8% of the initial envelope. A few, 
usually less than 10, SPH gas particles remain bound to the newly 
formed binary, been bound to either the newly formed WD, or the 
old WD. This explains why M\ can be slightly larger than M c \, 
and similarly why M 2 can exceed slightly 0.36 Mq. There is no cir¬ 
cumbinary envelope left in all simulations with the TEOS. In all 
our simulations, the final mass ratio q ranges between 0.88 — 0.90, 
consistent with the observational error for WD 1101+364. 

Energies. The total energy at the end of the simulation is dis¬ 
tributed in the “binding” energy of the gas bound to the binary, 
Fbountb the final orbital energy of the binary, £ 01 -b.fin, and the total 
energy of the unbound material at infinity, E“ unb : 

^-lot.fin Eorh.fin TEb oun d ^[oLunh• (it)) 

We have compared the initial and the final total energies and found 
that the error is less than 0.11% in all our simulations. 

£ tot,u„b is composed of E k ” n unb , E~ unb , E“ tunb - the kinetic, 
internal and potential energies of the unbound material, respec¬ 
tively. We note that £^ n unb is the dominant energy in the unbound 
material, though the internal energy of the unbound material at the 
end of the simulations is also non-negligible (see Table [2j. 

We present Eb OU nd f° r completeness, but the fate - accretion 
or ejection - of the several particles that remain bound to the bi¬ 
nary can not be resolved by the numerical method we use; on the 
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Table 2. Energies and masses 


Model 

^unb 

Mi 

m 2 

/7°° 
kin,unb 

poo 

^int,unb 

poo 

pot,unb 

poo 

^tot,unb 

^orb,fin 

Abound 

^tot,fin 

AEorb 

1.0RG0N 

0.663 

0.322 

0.360 

3.645 

0.473 

-0.025 

4.093 

-10.992 

-0.615 

-7.514 

-9.874 

I.0RG1N 

0.663 

0.322 

0.360 

4.123 

0.295 

-0.015 

4.403 

-11.278 

-0.582 

-7.457 

-10.220 

1.0RG2N 

0.663 

0.322 

0.360 

4.081 

0.543 

-0.024 

4.600 

-11.469 

-0.531 

-7.400 

-10.490 

1.2RG2N 

0.872 

0.323 

0.360 

4.604 

0.629 

-0.041 

5.192 

-15.504 

-0.639 

-10.951 

-14.159 

1.4RG2N 

1.079 

0.319 

0.360 

6.790 

0.907 

-0.094 

7.603 

-22.911 

-0.005 

-15.313 

-21.196 

1.5RG2N 

1.178 

0.319 

0.361 

6.089 

0.917 

-0.096 

6.910 

-25.484 

-0.026 

-18.600 

-23.469 

1.5RG2NP 

1.178 

0.320 

0.360 

7.415 

1.407 

-0.159 

8.663 

-26.969 

-0.366 

-18.665 

-24.958 

1.6RG2N 

1.274 

0.323 

0.362 

5.623 

1.812 

-0.440 

6.995 

-27.741 

-0.244 

-20.990 

-25.584 

1.6RG0S 

1.274 

0.323 

0.362 

5.603 

1.692 

-0.381 

6.914 

-27.228 

-0.309 

-20.623 

-24.987 

1.7RG2N 

1.370 

0.323 

0.366 

5.854 

2.042 

-0.417 

7.479 

-33.692 

-0.715 

-26.928 

-31.073 

1.7RG0S 

1.373 

0.323 

0.363 

5.032 

2.061 

-0.610 

6.483 

-32.417 

-0.466 

-26.400 

-29.723 

1.8RG2N 

1.478 

0.318 

0.362 

8.333 

1.675 

-0.371 

9.637 

-52.873 

-0.171 

-43.407 

-48.961 


M, nb , Mi, and Mi are the unbound, stripped RG core and old WD, in M @ . unb = £,-m“ nb v?/2, -E“ tunb = E,m" nb Ui, unb = Eand 
unb are kinetic, internal, potential and total energies, respectively, for the unbound material. 7s 0 rb,fin is the orbital energy after the CEE. E boun d is 
the total energy of the particles that remained bound to the binary. E tot g n is the total energy of all the particles. All energies are in I() 46 erg. 


timescale of our simulation they stay in an orbit within the RL of 
their stars. This energy includes the kinetic, internal, potential and 
recombination energies for these several SPH gas. 

We should clarify that E or b,nn does not have to match with 
the two-body approximation, namely £ or b = —GM\M2/(2a or \ i )- 
In the latter, the potential assumes a form 0 °c 1/r, while our 
code inc ludes the soften ed form as described in the appendix of 
iHernauist & Katz) (119891 ). When the separation between the two 
SPH special particles is more than two smoothing lengths, the po¬ 
tential reduces to the Keplerian form. However, this separation is 
less than two smoothing lengths for the point particles after the 
CEE, and the potential is softened accordingly. The difference in 
orbital energy between the two methods varies from about 3% (for 
1.0RG0N) to 14% (for 1.8RG2N), with the Keplerian values being 
closer to zero. The initial orbital energy, given by Equation[4] is the 
same as in the two-body approximation. 

Because energy is well conserved, we can equate Equations [6] 
and [TO] For that, we also use Equations |T] [5] and [7] If we neglect 
Eboundi we can re-write the conservation of energy in fractions of 
the change in the orbital energy: 

Obind + Okc + tCb- 1 - (H) 

We find that this is indeed the case in our simulations (see Table 3), 
and that the deviation from 1 is due to £bound : the maximum devi¬ 
ation occurs in 1.0RG0N (~ 7.8%) and the minimum in 1.4RG2N 
(~ 0.03%). Note that if Cfec = C£^ b = the previous equation re¬ 
duces to the standard energy formalism. However, values of both 
Orec and a“ b are non-negligible and comparable to a^ind- We em ‘ 
phasize that previously it had been anticipated only that abind is 
somewhat less than 1, and we provide new improved constraints. 
Unfortunately, this is not yet a final solution of the problem as cQ b 
can not be easily predicted for any system—a subject of our future 
studies. 

Orbital angular momenta. We find that the ejected material 
takes away more than 90% of the initial angular momentum of the 
binary. Values of 7 vary between 1.42 and 1.87. This large range of 
values unfortunately does not allow the obtained values of 7 to be 
useful for predicting the final parameters in a popul ation sy n thesi s 
for all possible DWD systems (for details, see llvanova et all2013h . 

Final orbital parameters. We find the final orbital separa¬ 
tion as a or b.fin = ( r a +r p )/2, where r p is the periastron, and r a is 



M [M 0 ] 


Figure 2. Final orbital periods versus initial mass of a RG. Triangles repre¬ 
sent non-synchronized RGs, and circles represent synchronized RGs. Note: 
to compare alike cases, we show only the r\ = 2 cases with our standard res¬ 
olution. Different rj’s give similar outcomes (see Table 3). 

the apastron. We ensure that these two quantities, r p and r a , do not 
change with time at the moment when we extract them from the 
simulations. We calculate the final orbital period Rorb.fin of the bi¬ 
nary from Kepler’s third law and the eccentricity of the post-CE 
orbit as e = (r a — r p )/(2a orb bn ). The latter is small in all the mod¬ 
els, showing that post-CE orbits are almost circular (in previous 
studies, where the ejection of the CE was incom plete, the fi nal ec¬ 
centricity was larger, 0.08 or more, e.g. lRicker & Taamll2012h . 

Figure [2] shows the final orbital periods plotted versus initial 
RG. We see that, as expected, the more massive the star is, the 
tighter the orbit gets. We find also that the final orbital period for the 
non-synchronized and synchronized cases are very similar (for the 
final state of the binary system, the only change due to syncroniza- 
tion was observed for the final eccentricity, albeit final eccentricity 
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Table 3. Orbital parameters 


Model 

•/orb, ini 

•/orb,fin 

7 

r p 

r a 

#orb,fin 

/orb, fin 

e 

^bind 

&rec 

<Cb 

1.0RG0N 

14.340 

1.188 

1.861 

2.015 

2.115 

2.065 

0.416 

0.024 

0.855 

-0.208 

0.431 

1.0RG1N 

14.741 

1.168 

1.868 

1.965 

2.074 

2.020 

0.403 

0.027 

0.827 

-0.201 

0.431 

1.0RG2N 

15.119 

1.157 

1.873 

1.947 

2.036 

2.000 

0.397 

0.022 

0.808 

-0.197 

0.440 

1.2RG2N 

16.262 

0.987 

1.670 

1.520 

1.532 

1.526 

0.264 

0.004 

0.871 

-0.192 

0.367 

1.4RG2N 

17.116 

0.759 

1.557 

1.070 

1.089 

1.080 

0.158 

0.009 

0.800 

-0.159 

0.359 

1.5RG2N 

17.062 

0.709 

1.512 

0.953 

1.003 

0.978 

0.134 

0.026 

0.879 

-0.172 

0.294 

1.5RG2NP 

17.062 

0.719 

1.511 

0.891 

0.924 

0.908 

0.122 

0.018 

0.815 

-0.148 

0.347 

1.6RG2N 

17.685 

0.678 

1.479 

0.880 

0.948 

0.914 

0.122 

0.037 

0.893 

-0.157 

0.273 

1ARGOS 

17.392 

0.690 

1.477 

0.912 

0.947 

0.930 

0.126 

0.019 

0.895 

-0.160 

0.277 

1.7RG2N 

17.151 

0.610 

1.449 

0.746 

0.771 

0.758 

0.092 

0.016 

0.922 

-0.140 

0.241 

1.7RG0S 

16.953 

0.624 

1.444 

0.776 

0.791 

0.784 

0.097 

0.009 

0.942 

-0.146 

0.218 

1.8RG2N 

14.932 

0.446 

1.417 

0.464 

0.493 

0.479 

0.047 

0.030 

0.902 

-0.096 

0.197 


The orbital angular momentum 7 0 rb,ini and ./orb,fin for the initial and final binary, respectively, in units of 10 52 g cm 2 s _1 . The parameter y is defined in Eq.[9] 
The closest and farthest orbital separations are r p and r a , respectively, while tf 0 rb,fin is the semimajor axis (all in Rq). The orbital period /orb, fin is given in 
days, and e is the eccentricity of the orbit. The energy fractions abind> °rec> and oQ b are defined in Eq.Q] Eq.[5] and Eq. [T| respectively. 


is small in all the cases). We conclude that our best progenitor for 
WD 1101+364 is a 1.4-1.5M 0 RG. 


4 CONCLUSIONS 

To understand the energy budget during a CEE leading to a DWD 
formation, we perform non- and synchronized 3D hydrodynamic 
simulations with two EOSs. We confirm that taking into account 
recombination energy leads to a full ejection of the RG’s envelope 
and the formation of a non-eccentric binary system, whilst if we do 
not take recombination energy into account, we obtain result simi¬ 
lar to previous studies and only half of the RG’s envelope is ejected. 
The most important consideration appears not to be the value of the 
available recombination energy, but where and when this energy is 
released. Indeed, ionized material forms the circumbinary enve¬ 
lope initially. Recombination then takes places there, while the cir¬ 
cumbinary envelope continues to expand. This results in the ejec¬ 
tion of the circumbinary envelope and effectively of all the com¬ 
mon envelope material. If instead the recombination energy had 
been released too early, the simulations would have ended up with 
unexpelled circumbinary envelope as in previous studies. In addi¬ 
tion, we find that considering a complete synchronization versus 
non-synchronized case does not change noticeably the final results. 

We introduce a modification of t he standard energy formalism 
dWebbinkl 1 984l ; [Livio & SokeJ 19881) , with the parameters describ¬ 
ing the use of the recombination energy and the unbound material 
energy. The first one can be found from initial stellar models, but 
the latter requires 3D simulations. For our set of models, a^ b has 
values from about 0.2 to about 0.44. However to generalize the re¬ 
sult and make it useful for population synthesis one needs to make 
a thorough parameter study; this is the subject of our future studies. 

As expected, we find that the more massive the parent RG star 
is, the tighter the final orbit gets. We do not find that the initial syn¬ 
chronization affects the final period but instead only changes the 
energy and angular momentum carried away by the ejecta, presum¬ 
ably shaping the post-CE nebula. We also find that our binaries end 
up with an eccentricity smaller than 0.04—a result that has been 
expected theoretically but not yet produced in simulations. 

We applied our met hod to the case of WD 1101+364, a well- 
known DWD (see iMarsiil 19951) . We inferred that its progenitor bi¬ 


nary could have been composed of a 1.4— 1.5M 0 RG and a0.36M Q 
WD companion, with Rorb.ini ~ 31 — 33 days. 
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